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CHAPTER 1 


INTRODUCTION 

The purpose of this thesis is to report the results of 
an empirical analysis of the duration of overshoots above an 
arbitrary level in a stationary gaussian stochastic process. 

The study was conducted for the Terrestrial Environment Branch, 
Aerospace Environment Division, Space Sciences Laboratory, 

George C. Marshall Space Flight Center, Alabama and was sup- 
ported under NASA contract no. NAS8-29286. The results 
presented in this thesis are applicable in the prediction of 
extreme properties of such natural processes as wind speed, 
ambient temperature, and sea state. 

1.1 Statement of the Problem 

The problem dealt with by this thesis concerns the 
frequency distribution of the duration of overshoots above an 
arbitrary level in a gaussian stochastic process with an 
exponential autocorrelation function. A gaussian stochastic 
process is a stochastic process which has a gaussian distri- 
bution at any point in time. Such a process is stationary if 
it has a mean independent of time and an autocorrelation 
function dependent only on the distance between successive 
time points. 

The problem is of general theoretical interest but 
numerical results are sparse. Rice (1945) did the fundamental 
work in studying both the crossing frequency and duration. 



Favreau, Low, and Pfeffer (1956) studied the duration problem 
and hypothesized the negative exponential density for duration 
time--this hypothesis was disproved by Longuet-Higgins (1962). 
The recent texts by Kramer and Leadbetter (1967) and Kuznetsov 
(1965) give excellent summaries of the work in this area. To 
the author's knowledge this is the first investigation of 
duration times that is based on an extensive simulated data 
set. 


1.2 Organization of the Analysis 

Chapter 2 contains a description of the development of 
the simulation model and its underlying assumptions. In 
Chapter 3, the analysis of the model and our conclusions are 
presented. Chapter 4 contains a discussion of the application 
of the results of this analysis to specific problems con- 
cerning atmospheric variables. 
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CHAPTER 2 


. ? 

Simulation Model 

Our first objective was to de.'ign a mathematical model 
of a stationary gaussian stochastic process with an exponential 
autocorrelation function- To do this, we made the following 
assumptions : 

1) The process had a multivariate normal distribution. 

2) The autocorrelation function, R(t), was exponential, 
i .e . , R ( t) = exp (-R | T | ) . 

3) The process was stationary, i.e., R(t^,tj) = R ( t) , 
where t = 1 1^ - t ^ | . 

4) The expected value of a random variable X at time t 
was zero, i.e., E[X(t)] = 0. 

5) The variance covariance matrix l was symmetrical and 
positive definite. 

In the remainder of this chapter, we shall use X(t) to 
denote a stochastic process satisfying the above conditions. 

The process was considered over a time interval of 
[ 0 , 991 , each realization consisting cf 100 equally spaced 
points in the interval. In the case of a specific application, 
the interval of interest may be [0,T]. X(t) would then be 
sampled at 100 equally spaced points in the interval. These 
points would be t Q , , . . . , t 99 , where t i = (f^)T, i = 0,1, ...99. 

The autocorrelation parameter R would also require modification. 

Odell (1971) presents the method of simulation used here. 

A summary of that technique follows. 
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Let X denote the vector 

[X(t Q ) , 

x(t 1 ) 

, . . . ,X(t 

W r ' 

where ' denotes matrix transposition. 

Let 

Z » (C..) » E(X-X'), 

where E is the expectation operator. 

E is 

, of course, the 

variance covariance 

matrix. 

Thus for 

each 

element 

of 

E, o i; ., 0 < 

j I 

99 , we have 




Oi j = E[X(t i ) • X ( t 

j)l - R(t 

.,t.) * R(t) , ' 

where t 

* iv s'- 

Hence we see that 
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R(0) 



By assumption, X satisfies a multivariate normal distribution 
with mean y = 0 and covariance matrix Z. We write X 'v N(y,E). 
From Odell (1971, pg. 37) we have the following: 

Theorem: If the 100 x 1 vector Y ^ N(y,E), 

| is a fixed 100 x 1 vector, and 

V = AY + then 

V ^ N(Au + X» A 1 A '> • 

We generated the vector Y ^ N(0, I), where I is the identity 
matrix. We factored the covariance matrix E into AA' by the 
Crout method. Thus, if V * AY, it follows that V ^ N(0,E). 

The vector V constituted a realization of the process X(t). 

The vector Y * (y 0 »y^» • • • rYgg) was formed by generating a 
sequence of 100 independent standard normal variates. 

Hamming (1962) provides the technique for generating Y. This 
method is summarized below. 
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An approximation to normally distributed random numbers 
may be obtained from a sequence x k , k = 1, 2,...,K of uni- 
formly distributed random numbers by the formula 
K 

E x. - K 

y. = k = 1 * . To utilize this method, the value of K 

i — 

/ K/12 

12 

was fixed at 12, reducing the formula to y . - I x, - 6. 

1 k= J K 

Thus, we were able to produce a sequence 

Y = (y^ ,y i , . . . ,y^g ) of standard normal variates with mean 0 
and unit variance. Hence, we had Y ">■ N(0,T). 

For each autocorrelation function, we qenerated 250 
realizations which required 250 random vectors Y_j , i = l , . . . , 250 , 
each of which required 100 standard normal variates. Each 
normal variate required 12 uniformly distributed random num- 
bers. Hence, for each autocorrelation function, 

250 x 100 x 12 = 300,000 uniformly distributed random numbers 
were required. The power residual method was used to generate 
these random numbers, each one being normalized to (0,1). 

That is, 

r = normalized (s ), where s = X*r ,. 
n n n n- 1 

The computer used for the generation was a Univac Series 

31 

70/46 which has an integer capacity of 2 - 1. 

In order to transform each vector Y^ into a realization 
, via A Y^, it was necessary to factor E into A • A ' . Since 

E is a symmetrical, positive definite matrix, it can be 
factored into a lower triangular matrix and its transpose. 
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The factorization was accomplished using the Crout factorization 
method, Odell (1971, pg . 38). This method is summarized in 
Carter and Madison (1973) . 

The autocorrelation function R(t) = exp(-B|T|) determines 
the degree of association between pairs of values of X(t). We 
simulated the process for twelve autocorrelation functions by 
assigning 8, the autocorrelation parameter, the values 
.005, .010, .025, .050, .075, .10, .25, .50, .75, 1.0, 2.5, 5.0. 
With 8 = .005, the process had a correlation above .61 through- 
out and .99 for successive pairs of points. For 8 = 5.0, 

R(2) = .00004539, meaning that the values were essentially 
independent after two time intervals. The range of 8 values 
used proved to be appropriate for analysis purposes. 
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CHAPTER 3 


Data Analysis 

With the data sets generated, as described in Chapter 2, 

the first problem was that of determining the durations of 

overshoots above selected levels. From Sveshnikov (1966), we 

see that in a continuous random process X(t), X:D*R, two co r - 

ditions must be met in order for an overshoot above a level a 

to occur in the interval (t,t + dt) . They are that at a given 

time t we have X(t) ' a, while X(t + dt) > a. We say that an 

overshoot occurs between t^ and t provided X(t^) a, 

X(tj) ' a, and X ( t ) N a for all tr (t^ , t ^ )H d . Tn the discrete 

capo, we see that, the two conditions required for an overshoot 

to occur above a level A are that X(t n ) ^ A, while X(t + j) ' A, for 

some n. We say that an overshoot occurs between t and t 

n m 

provided X(t ) ' A, X(t ) ' A, and X(t.) A for n ' i ' m. 
n — m — l 

Figures 1 and 2 illustrate an overshoot in the continuous case 
and one in the discrete case, respectively. 



Tn trie continuous case, the duration of the overshoot between 

t^ and tj would clearly be (t^ - t^) . In the discrete case, 
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to determine the duration we can only count the number of points 

t. between t and t such that X(t ) > A, while X(t ) < A, and 
i n m i n — 

X(t ) < A, which is m - n - 1. 
m — 

Eac’' realization in our data set consisted of 100 points 

X(t Q ), X(t X(t<^). So for a given level A, we located 

each overshoot, that is we found i such that X(t. ) f h, while 

l-l ~ 

X(t J > A, and then counted the points until X(t.) <_ A; the 

duration was j-i. (0 < i < j 99). The question soon 

surfaced of whether to count as overshoots those points that 

began above A, X(t^) > A, or which ended above A, X(t^) > A. 

Noting that for small 8 values the system has a let of memory, 

if the process began with X(t^) > A, it was probable that it 

would stay above A for a long time. Similarly, if the process 

exceeded A somewhere between t^ and t , it was likely to stay 

0 99 

bove A up to and including t . For larger 8 values, we did 
not feel that counting the'e so called TO and T99 overshoots 
would appreciably affect our results. We therefore decided to 
include these </.ypes of overshoots in our counting. 

We counted overshoots above levels A = 0, .5, .75, 1, 1.25, 
1.5, 1.75, 2. 2.5, 3, 3.5, 4. Recalling that the process has 
a mean 0^ and unit variance, we were equivalently considering 
overshoots above the mean, .5 standard deviations above the 
mean, .75 standard deviations above the mean, etc. In applications, 
for a process with mean u and variance a 2 , the overshoots above 
level A would be equivalent to overshoots above the value of 
A • o + y . 
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The value 3.0 was selected as the upper limit of the A 
levels for analysis purposes. In the independent case, only 
.13% of the values would be above 3.0, and in the more cor- 
related cases, the number of points above 3.0 and thus, the 
number of overshoots, would likely decrease. Zero was selected 
as the lower limit. 

The durations of the overshoots for each pair, level A 
and autocorrelation parameter 6, were determined and the sample 
mean X and standard deviation were computed as usual. 

Table 1 contains the means for each pair (A,B) , and Table 2, 
the standard deviations. 
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Our next step was to plot the mean X vs. the crossing 
level A and the standard deviation S vs. the crossing level. 
Several of these curves are shown in the figures 3 and 4. 

Our search for estimating equations thus began. The curves 
indicated an exponential trend, becoming flatter and con- 
verging as the 6 value increased. So we graphed the same 
data on log paper. These plots were much strai hter, but 
there was still some curvature left in them. Suspecting an 
exponential characteristic again, we tried log-log paper; the 
results were not satisfactory. We decided that the curvature 
left from the log paper could probably be taken care of with 
a quadratic function. 

We ran polynomial regression analyses, BMD (1973) , for X 
as a function of A and lnX as a function of A. We concluded 
that we could best fit the curves using lnX as a quadratic 
function of A. We did the same for the standard deviation S 
and concluded that our best fit would be obtained by con- 
sidering S as a quadratic function of A. We thus had the 
following: _ 

lnX = a^ + a^A + ajA 2 and (III.l) 

S=b + b A + b A 1 , where the coefficients 
0 1 2 

Sf jbf , i=0 , 1 ,2 were dependent on 6, the 
autocorrelation parameter. 

Table 3 gives the simple correlation coefficients for 
equations (III.l) for each 8 value. Of course, the presence 


12 



of the quadratic term did, in every case, increase the 
correlation between observed and predicted values. 
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TABLE 3 

SIMPLE CORRELATION COEFFICIENTS 


Beta 

For Means 

For Standard 
Deviations 

.005 

-n 95346 

-0.96940 

.010 

-0.90434 

-0.94986 

.025 

-0.98143 

-0.98375 

.050 

-0.81609 

-0.85630 

.075 

-0.99292 

-0.95819 

.100 

-0.99394 

-0.95653 

.25 

-0.97928 

-0.91802 

.50 

-0.93697 

-0.91145 

.75 

-0.95889 

-0.94079 

1.0 

-0.94336 

-0.94624 

2.5 

-0.87894 

-0.91315 

5.0 

-0.87083 

-0.91963 
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It is clear from the table that the models (equations III.1) 
are adequate for each 8 value. 

At this point, for each B_, we had the coefficients 
aQ,a lf a 2 for the mean and bQjb^b^ for the standard deviation. 

To determine the relationship between P and the constant terms 
for the men' i, we plotted ag as a function of 6. There was 
little doubt that this curve behaved in an exponential manner. 
The same phenomenon was observed in plots of other coefficients 
as functions of R for both the mean and standard deviation. 

We considered the coefficients as functions of InR, and then 
ran polynomial regression analyses of degree two and three. 
Choosing the third degree, we obtained very good fits. Hence 
we have the following: 

2 3 

a^ = Cq + c^(lnR) + C 2 (lnR) + c^(lnp) , and 

o , (III. 2) 

b^ = d Q + d^ ( In R) + d^lnp; + d^Qnp) , 

for i = 0,1,2. 

The coefficients c^ and d^ , i**0, 1,2,3, are presented in Table 

4 for each coefficient a.,b., i»0,l,2. 

i l 
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b Q 1.88589 -1.02364 0.56796 -0.09535 

b^ -1.17242 0.74806 -0.38702 -0.01600 

b 0.20194 -0.16723 0.07649 0.01833 

2 
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Using the last set of coefficients, we were able to 
predict the values a^ , a^ ,b^ ,b^ ,b^ for each B value. Hence, 
for each B value and each A level, using (III. 2) and substituting 
into (III.l), we could predict X and S as functions of P and A. 

In summary, we have: 

_ 2 

X predicted = exp[a^(lnP) + a^(ln8)*A + a 2 (lnp)*A ], and 

S predicted = b^dnB) + b^(ln8)*A + bjdnBJ'A^, where 
a^dnS) is a function of ln8, etc. 
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Figures 3 and 4 picture some of the observed curves 
together with the corresponding predicted curves for the 
means and standard deviations, respectively. As can be seen 
from the plots, for both the means and the standard deviations 
the predicted curves match the observed curves quite well. 

For the means, the predicted values at A=0 are, consistently, 
slightly less than the observed values. Throughout the range 
of A values, the predicted curves for the means and standard 
deviations appear to fit equally well. Overall, the fits were 
much better than one would normally expect in such an analysis 
An interesting phenomenon was observed while viewing 
the plots of the linear and quadratic coefficient curves. 
Figures 5 and 6 picture the plots of the linear coefficients 
a^ vs. the quadratic coefficients a 2 for the means and stan- 
dard deviations, respectively. There is definitely a high 
correlation between these values. The author has no concrete 
explanation to proffer. It may be simply coincidence, it may 
be related to the polynomial regression program used, or it 
may be a common occurence in such a situation. From a 
sl-atistical standpoint, the similarity means that given the 
linear coefficient, we are also given prior knowledge of the 
quadratic coefficient. Thus the linear and quadratic terms 
would not be totally independent. At any rate, we do not 
believe that our results were significantly affected by this 
event . 
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W<~> constructed frequency tables for each A level wit.hin 
each B, and Hien made plots of observed cumulative frequencies. 
Some of these are presented in Appendix A. The cumulative 
plots resembled some type of gamma distribution. We decided to 
first try the exponential distribution which is a special case 
of the qamma. We used as estimates of theta 1/observed mean, 
1/observed standard deviation, 1/predicted mean, and 1/predicted 
standard deviation. The results were not altogether bad, but 
the estimated frequencies were consistently high for low A 
levels. We next tried the gamma distribution whose density 
function is 

b a 

f(x;a,b) ~ rjaT * x * ^ * e > f°t* X'O, 0 elsewhere. 

The qamma distribution is a two parameter distribution (the 
exponential is a special case) and therefore provides us a 
richer family of curves to use. The gamma distribution, its 
many uses and properties, is well described in Johnson and 
Kotz (1970). The mean u and variance <-|2 are given by: 

u - a/b and o* = a/b . Noting that * b and Ji ? * a, when 

2 

given the process mean X and variance S , we can, by the 

method of moments, easi ly estimate the parameters a and b for 

the gamma distribution pdf tor any of our processes. 

Using the corrected mean (X - 1/2) , we utilized the gamma 

— 2 2 — 2 

distribution with parameters a ■ (X -.5) /S and b * (X -.5)/S 
to find the predicted probabilities of overshoots of duration 
1,2,3, etc. For example, for a particular B value and A level, 
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we used either the predicted or observed . 1 .ies of the mean 
and vari.jr 1 '* to calculate a and b. Then. *- determine the 


probabi 1 i t y of the occurence of an overshoe’ of duration 1, 
we integrated as follows: 


Probability 


; 3/2 j£_ 

i n r(a) 



• dx . 




(x 



We then used the Kolmogorov-Smi rnov goodness of fit test to 
determine our degree of accuracy. The results were quite 
favorable. Appendix B presents the results of the distribution 
fitting for the A and 6 values presented in Appendix A. 

Table 5 gives an overall scenario of the goodness of fit 
results. Upon viewing the table several points become clear. 
The gamma distribution gives excellent results for most A and B 
combinations. Especially important is the consistently excel- 
lent results for A values above 1.5--the range of A values 
utilized most frequently. The combination of A and B values 
where the gamma model fails is quite interesting. v.hen 
comparing the predicted mean and variance with the observed 
values one notes little difference and, in fact, the gamma 
model using the observed values was not satisfactory either. 

There are three, possibly related, reasons why the model 
fails in these cases. The first is the very conservative 
nature of the Kolmogorov-Smirnov test when the model parameters 
are estimated from the data (indirectly using the predicting 
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equations) and the data are discrete although the model is 
a continuous one. The second reason is the possible inadequacy 
of tie qamma model. A survey of Appendix B shows the rejection 
always occur.- at the first frequency because the qamma model 
underestimates the frequency. This suggests that a more com- 
plex model (perhaps a four parameter gamma) might be more satis- 
factory. However, the third reason is, in the author's opinion, 
the real culprit. Quite simply, the sample sizes are too 
large. Tt is an interesting phenomenon in statistics that in 
ooodnoss of fit problems too much data poses as many problems 
as does too little data. Lancaster (1969, pp. 174-75) discusses 
this problem and recommends utilization of moderate sample sizes. 
In no case where the fit is bad would it have been so if the 
same cumulative proportions had been observed with sample sizes 
under 500. 

Of course, to actually varify that this phenomenon 
actually causes the problem would be difficult, if not impos- 
sible. However, it is pertinent to note that the difficulty 
starts when the sample sizes become large while the mean and 
variance predictions remain quite accurate and, as one can 
observe in Appendix A, the actual shape of the sample cumu- 
lative distribution curves remain consistent. 
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TABLE 5 

Condensed Summary of Fitted Models 
Using Predictive Equations 

*Good Fit at a < .01 yt Bad Fit due to mean overestimation 

**Gocd Fit at a < .05 yt Bad Fit due to mean underestimation 

0 No data a 2 +Bad Fit due to variance overestimation 

o 2 *Bad Fit due to variance underestimation 
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CHAPTER 4 

Appl i.cations 

This chapter presents the procedures to utilize the 
results obtained in Chapter 3. The computer program developed 
to do the necessary calculations is presented in Appendix C. 

The analysis was conducted using a standard normal process 
2 

(m = 0, o =1) and duration times above a level A were counted. 
However, a simple transformation permits the use of any mean 
and variance, and duration times below a specified level can 
be addressed using the symmetry of the normal distribution. 

A final note involves the time intervals. In most 
processes the sampling and estimation procedures use data 
gathered at regular intervals and the results are, to some 
extent, dependent on the interval widths. Questions must be 
posed in this context, e.g., it is not reasonable to request 
the probability of a duration time exceeding five seconds when 
the autocorrelation parameter is estimated based on samples at 
one minute intervals — the autocorrelation parameter estimate 
must also be based on (or modified to correspond to) data 
gathered at one second intervals. 

Our simulation process, hence the 8 values, used unit 
intervals corresponding to serial correlation coefficients of 
lag 1, lag 2, etc. Should a B value using the actual time 
intervals be given, it must be modified to represent the auto- 
correlation parameter of a process sampled at unit intervals 
and the time period of these intervals must agree with the 
time period referred to in the problem under consideration. 
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This transformation procedure is summarized as follows. Assume 
a value B' has been calculated using intervals t' = 0 , L , 2L , 3L . . . 
giving R(x') = exp(-6'jx|). This corresponds to the function 
R(x) = exp (-6 "h | -t| ) where t= 0,1,2,... . Thus 8 = p "h would 
be the autocorrelation parameter estimate corresponding to 
intervals of unit length. 

As an example, assume a wind speed process where 
y = 24 m/s, o = 8 m/s, and 6' = .026 based on samples measured 
at 10 minute intervals. We desire to calculate the probability 
of exceeding 30 m/s for 3.2 minutes. The 6' value is for 
t" = 0, 10 sec., 20 sec.,... . The quantity 

8 = 6 'h = .026 * 10 = .26 is for x = 0,1 sec., 2 sec.,... and 
is the value we would utilize in calculating the probability. 

The computer program required the following inputs: 

1) Process mean y, 

2) Process standard deviation o, 

3) Process autocorrelation parameter (modified 

to reflect intervals of unit lengtl j.f necessary) , 

4) The crossing level L under consideration, and 

5) The number of time units x (perhaps fractional) 
that the duration lasts. 

The output gives 

1) The y, o, 8, and L values, 

2) The adjusted crossing level A (= |L-y|/o), 

3) Predicted mean and variance of the process, and 

4) Pr (Exceeding L for x time units). 
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Referring to the example in the preceding paragraph the 
inputs would be: 

1) v = 24, 

2) o=8, 

3) 0 = .26, (the .026 value was modified), 

4) L = 30, and 

5) x = 3.2. 

Table 6 gives the computer results corresponding to these 
inputs . 

The question concerning the applicability of these results 
to nonnormal processes remains. Carter and Madison (1973, pp. 
35-36) discussed this problem at length. Basically the procedure 
outlined above will give satisfactory results when the parent 
population is normal appearing - the case with many nonnormal 
processes for certain ranges of values for the population 
parameters. The estimating procedure developed in Chapter 3 
utilizes a unique property of the normal process - independence 
of the sample means and variances. Likely the nonnormal process 
must be of a form to make this assumption plausible for our pro- 
cedure to work satisfactorily. 
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TABLE 6 

PREDICTION PROGRAM OUTPUT 
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APPENDIX A 

SAMPLE CUMULATIVE DISTRIBUTION CURVES 
FOR AUTOCORRELATION PARAMETERS: 

.025, .075, .25, .75, 1.0, 2.5, 5.0 
AND A-LEVELS: 0., .75, 1.25, 1.75, 2.5 
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APPENDIX B 

KOLMOGOROV-SMIRNOV GOODNESS OF FIT TESTS 
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KOIMOCOROV-SMIRNQV GOODNESS 0 F f IT TEST 
using predicted values 


A-LEVEL • 0.000 AUTOCORRELATION PARAMETER • 0.025 

OBSERVED MEAN ■ 12.15896 OBSERVED VARIANCE ■ 5*2,92211 

PREDICTED MEAN ■ 11.605*1 PREDICTED VARIANCE • 530, *0380 


GAMMA OIST 

PARAMETERS! 

ALPHA ■ 

0,37327 BETA 

• 0,03361 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

OBSERVED 

OBSERVED 

expected 

DIFFERENCE 

1 

306 

0.28281 

0.31*13 

*0,03132 

2 

1*0 

0. *1220 

0,*032* 

0,00896 

3 

72 

0. *787* 

0.*6*98 

0.01377 

* 

5* 

0.52865 

0.5131* 

0,01551 

5 

*6 

0.57116 

0.55287 

0,01829 

6 

36 

0.60*** 

0.5*672 

0,01771 

7 

*2 

0.64325 

0.61619 

0,02707 

8 

28 

0.66913 

0.6*223 

0,02690 

9 

15 

0.68299 

0.66550 

0,017*9 

10 

19 

0*70055 

0.686*9 

0,01*07 

11 

11 

0.71072 

0,7055* 

0,00518 

12 

23 

0.73198 

0.77295 

0,00903 

13 

12 

0.7*307 

0.73892 

0,00*15 

1* 

9 

0.75139 

0,7536* 

-0,00225 

15 

1* 

0.76*33 

0.7672* 

-0,00292 

16 

8 

0.77172 

0.77986 

*0 , 008 1 * 

17 

12 

0.78281 

0.79159 

*0,00878 

18 

10 

0.79205 

0.80252 

-0,010*7 

19 

15 

0.80591 

0.81273 

-0,00682 

20 

10 

0.8iSl6 

0.82228 

*0,00712 

OVER >0 

200 

1,00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 1082 
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KflLMOGORDV-SMlRNOV GOODNESS OF FIT TEST 
USING PREDICTED VALUES 


A-LEVEL ■ 0.750 AUTOCORRELATION PARAMETER • 0.025 

OBSERVED MEAN » 7.79122 OBSERVED VARIANCE • 1*7,37265 

PREDICTED MEAN • 8.12091 PREDICTED VARIANCE ■ 153,98528 


GAMMA DIST 

PARAMETERS! 

ALPHA ■ 

0.37717 BFTA 

• 0,04949 


number 

cumulative 

CUMULATIVE 


DURATION 

observed 

OBSERVED 

expected 

difference 

1 

239 

0.31782 

0,3*727 

-0,03945 

2 

114 

0.469*1 

0.45791 

0,01150 

3 

62 

0.55186 

0,52666 

0,02520 

4 

4* 

0.6i037 

0.57952 

0,03085 

5 

37 

0.65957 

0.62248 

0,03710 

6 

29 

0.69814 

0.65853 

0,03961 

7 

23 

0.72872 

0.60943 

0,03929 

8 

18 

0.75266 

0.71632 

0,03633 

9 

9 

0.76463 

0.74000 

0,02463 

10 

14 

0.78324 

0.76101 

0,02223 

11 

14 

0.80186 

0.77981 

0.02205 

12 

10 

0.81516 

0,79670 

0,01845 

13 

12 

0.83112 

0.81197 

0,01914 

14 

12 

0.84707 

0.82582 

0,02125 

15 

10 

0.86037 

0,83843 

0,02194 

16 

8 

0.87101 

0,84994 

0,02107 

17 

6 

0.87899 

0.86047 

0,01852 

18 

3 

0.88298 

0.87014 

0,01284 

19 

* 

0.88830 

0.87902 

0,00928 

20 

3 

0.89229 

0,88720 

0,00508 

OVER 20 

81 

1.00000 

1.00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 752 
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KfllMOGOROV-SMlRNOV GOODNESS 0* FIT TEST 

USINO predicted values 


A-LEVEL ■ 1.290 AUTOCORRELATION PARAMETER ■ 0.025 

OBSERVED MEAN • 6*21862 OBSERVED VARIANCE • 95,97799 

PREDICTED MEAN • 6.32)18 PREDICTED VARIANCE • 63,7*119 


GAMMA OIST 

PARAMETERS! 

ALPHA • 

0.40493 BETA 

• 0,06954 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

observed 

OBSERVEC 

EXPECTED 

DIFFERENCE 

1 

164 

0.38588 

0.37552 

0,01036 

2 

61 

0.52941 

0.48763 

0,04178 

3 

38 

0.61882 

0.56381 

0,05501 

4 

22 

0.67059 

0.62177 

0,04882 

5 

17 

0.71059 

0.66825 

0,04234 

6 

15 

0.74588 

0.70670 

0,03918 

7 

9 

0.76706 

0.73916 

0,02790 

8 

10 

0.79059 

0.76696 

0,02363 

9 

5 

0.80235 

0.79102 

0,01133 

10 

13 

0.83294 

0.81203 

0,02091 

11 

6 

0.84706 

0,83049 

0,01657 

12 

7 

0.86353 

0,64680 

0,01673 

13 

7 

0.88000 

0.86128 

0,01872 

14 

4 

0.88941 

0.67418 

0,01523 

15 

4 

0.89882 

0.88571 

0,01311 

16 

0 

0.898B2 

0.89605 

0,00277 

17 

2 

0.90353 

0.90534 

•0,00181 

18 

2 

0.90824 

0.91371 

•0,00547 

19 

2 

0.91294 

0.92126 

•0,00832 

20 

2 

0,91765 

0.92809 

•0,01044 

OVER 20 

35 

1.00000 

i; 00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 629 
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KniMOGOROV-SMIRNQV GOODNESS QF FIT TEST 
USING PREDICTED VAlUFS 


A-LEVEl ■ 1.750 AUTOCHRRElAT ION PARAMETER • 0.025 


OBSERVED MEAN * 5.10377 OBSERVED VARIANCE ■ 37,66846 

PREDICTED MEAN « 4.87360 PREDICTED VARIANCE ■ 40,40451 


gamma dist 

PARAMETERS! 

ALPHA • 

0.47292 BETA 

« 0,10808 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

observeo 

observed 

fxpectfd 

DIFFERENCE 

l 

76 

0.35849 

0.38100 

-0,02251 

2 

35 

0.52358 

0.51154 

0,01205 

3 

19 

0.61321 

0.60003 

0,01318 

4 

12 

0.66981 

0.66632 

0,00349 

5 

7 

0.70283 

0,71830 

-0,01555 

6 

11 

0.75472 

0,76037 

-0,00566 

7 

9 

0.79717 

0,79488 

0,00229 

8 

2 

0.80660 

0.82359 

-0,01698 

9 

3 

0.82075 

0,84771 

-0,02695 

10 

1 

0.02547 

0.868J2 

-0,04265 

U 

6 

0.85377 

0,88550 

-0,03172 

12 

7 

0.88679 

0,90036 

-0,01357 

13 

3 

0.90094 

0.91313 

-0,01219 

14 

2 

0.91038 

0.92414 

-0,01376 

15 

1 

0.91509 

0.93365 

-0,01855 

16 

3 

0.92925 

0,94189 

-0,01264 

17 

2 

0.93868 

0.94905 

-0,01037 

18 

0 

0.93868 

0,95527 

-0,01659 

19 

3 

0.93283 

0,96070 

*0,00787 

20 

0 

0.93283 

0,96344 

•0,01261 

OVER 20 

10 

1.00000 

1,00000 

0,00000 

total number 

OF OVERSHOOTS ■ 212 
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KOIHOGORUV-SMIRNOV GOODNESS OF FIT TEST 

using predicted values 


A-LEVEL ■ 2.500 AUTOCORRELATION PARAMETER ■ 0.025 


OBSERVED MEAN - 2,6*286 OBSERVED VARIANCE - 6,9*988 

PREDICTED MEAN • 9.2*129 PREDICTED VARIANCE • 9,97086 


gamma dist 

PARAMETERS! 

ALPHA ■ 

0.80192 BETA 

• 0,29253 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

observed 

observed 

f),PECTFD 

DIFFERENCE 

1 

22 

0.52381 

0.38289 

0,17092 

2 

6 

0. 6h667 

0,5*617 

0 i 120*9 

3 

5 

0.78571 

0,67599 

0,10972 

* 

1 

0.80952 

0.76651 

0,0*301 

5 

2 

0,8571* 

0.83075 

0,02639 

6 

1 

0.88095 

0.876*1 

0,00*1* 

7 

* 

0.97619 

0,91007 

0,06612 

R 

0 

0.97619 

0.98*19 

0,0*200 

9 

0 

0.97619 

0.9*175 

0,02*** 

10 

0 

0.97619 

0.96*57 

0,01162 

11 

0 

0.97619 

0,97395 

0,0022* 

12 

0 

0,97619 

0.9*083 

•0,00*6* 

13 

0 

0.97619 

0.9*587 

•0,00968 

1* 

1 

) .00000 

0.9*958 

0,010*2 

15 

0 

1.00000 

0.99231 

0,00769 

16 

0 

1.00000 

0.99*33 

0,00567 

17 

0 

1.00000 

0.995*1 

0,00*19 

18 

0 

1.00000 

0.99690 

0,00310 

19 

0 

1.00000 

0,99771 

0,00229 

20 

0 

1.00000 

0.99830 

0,00170 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • *1 
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KfUMOGOROV-SMlRNOV GOODNESS OR RIT TEST 
USING PREDICTED VAlUFS 


A-LEVEL ■ 0.000 AUTOCORRELATION PARAMETER • 0,075 


OBSERVED MEAN ■ 7,73678 OBSERVED VARIANCE • 109,75586 

PREDICTED MFAN • 7,32331 PREDICTED VARIANCE ■ 100,10458 


gamma dist 

parameters l 

ALPHA ■ 

0.46509 BETA ■ 

0,06816 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVED 

OBSERVED 

expected 

oifferenci 

1 

463 

0.28811 

0.31691 

-0,02880 

2 

231 

0.43186 

0 , 4 2 8*6 

0,00350 

3 

126 

0.51027 

0,50670 

0,00357 

4 

101 

0.57312 

0,56763 

0,00549 

5 

70 

0.61668 

0,61732 

—0 , 00064 

6 

67 

0.65837 

0,66899 

-0,00062 

7 

62 

0.69695 

0.69456 

0,00239 

R 

52 

0.72931 

0,73574 

0,00397 

9 

42 

0.75544 

0,75232 

0,00322 

10 

29 

0.77349 

0,77588 

-0,00239 

11 

24 

0.78843 

0,79682 

-0,00840 

12 

27 

0,80523 

0,81546 

-0,01023 

13 

20 

0.81767 

0.83210 

-0,01443 

14 

21 

0.83074 

0,84702 

-0,01628 

15 

12 

0.83821 

0,86044 

-0,02223 

16 

19 

0.85003 

0.87253 

-0,02250 

17 

23 

0.86434 

0.86345 

-0,01911 

18 

18 

0.87354 

0.89334 

*0,01779 

19 

9 

0.88114 

0.90230 

-0,02115 

20 

25 

0.89670 

0 , 91044 

-0,01374 

OVER 20 

166 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 1607 
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KOLMOGOROV- SMIRNOV GOODNESS 0? FIT TEST 
USING PREDICTED VALUES 


A-LEVEL ■ 0.750 AUTOCORRELATION PARAMETER • 

OBSERVED MEAN • *,65253 OBSERVED VARIANCE ■ 

PREDICTED MEAN ■ 4,8253* PREDICTED VARIANCE • 

GAMMA DIST PARAMETERS! ALPHA ■ 0.*3669 BETA • 



NUMBER 

CUMULATIVE 

CUMULATIVE 

DURATION 

OBSERVEO 

observed 

fXPECTFD 

1 

*30 

0.357** 

0.391*5 

2 

200 

0.52369 

0.520*8 

3 

127 

0.62926 

0,607*0 

* 

80 

0.69576 

0.67233 

5 

67 

0.751*5 

0.72322 

6 

5* 

0.7963* 

0.76*25 

7 

33 

0.82377 

0,79795 

8 

33 

0,85121 

0,82600 

9 

25 

0.87199 

0.8*957 

10 

18 

0,88695 

0,86953 

11 

19 

0.9027* 

0.8865* 

12 

9 

0.91022 

0.9O1U 

13 

1* 

0.92186 

0,91363 

l* 

7 

0.92768 

0.92*** 

15 

9 

0.93516 

0,93379 

16 

10 

0.943*7 

0.9*191 

17 

6 

0,948*6 

0,9*897 

18 

10 

0.95677 

0.95512 

19 

5 

0.96093 

0,960*9 

20 

7 

0,96675 

0,96519 

OVER 20 

*0 

1,00000 

1,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 1209 


0,075 

35,19763 

*0,96538 

0,10559 


DIFFERENCE 

- 0 , 03*01 

0,00321 

0,02186 

0 , 023*3 

0,02823 

0,03209 

0,02582 

0,02521 

0 , 022*2 

0 , 017*2 

0,01620 

0,00912 

0,00823 

0 , 0032 * 

0,00137 

0,00156 

- 0,00051 

0,00165 

0 , 000 ** 

0,00156 

0,00000 
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KOLMQGOROV-SMlRNQV GOODNESS op pit test 
USING PREDICTED VALUES 


A-LEVEL ■ 1,250 AUTOCORRELATION PARAMETER • 0,075 


OBSERVED MEAN • 3.58438 OBSERVED VARIANCE ■ 19,36970 

PREDICTED MEAN • 3,71363 PREDICTED VARIANCE • 20,39255 


gamma dist 

PARAMETERS! 

ALPHA m 

0.50706 BETA 

• 0,15769 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

observed 

OBSERVED 

expected 

DIFFERENCE 

1 

304 

0.42399 

0.41976 

0,00423 

2 

123 

0.59554 

0,56756 

0,02798 

3 

79 

0.70572 

0,66461 

0,04111 

4 

53 

0.77964 

0.73461 

0,04503 

5 

25 

0.81430 

0,70736 

0,02715 

6 

30 

0.85635 

0 . 8 2 8 i 3 

0,02821 

7 

22 

0.88703 

0,86019 

0,02684 

8 

12 

0.90377 

0.80570 

0,01806 

9 

13 

0.92190 

0.90619 

0,01571 

10 

6 

0.91026 

0,92274 

0,00752 

11 

8 

0.94142 

0,93620 

0,00522 

12 

6 

0.94979 

0,94719 

0,00260 

13 

1 

0.93119 

0.95620 

-0,00501 

u 

8 

0,96234 

0,96361 

-0,00126 

15 

4 

0.96792 

0,96971 

-0,00179 

16 

3 

0.97211 

0.97476 

-0,00266 

17 

2 

0,97490 

0,97894 

• 0 , 00405 

18 

3 

0.97908 

0.90241 

•0,00333 

19 

1 

0.98047 

0,90529 

•0,00482 

20 

0 

0.98047 

0,90749 

•0,00722 

OVER 20 

14 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOT* • 717 
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KOLMOGOAOV-SMlRNOV GOODNESS OF fit test 
USINO PREDICTED VALUES 


A-LEVEL ■ 1.750 AUTOCORRELATION PARAMETER • 0.075 

OBSERVED MEAN « 5,00961 OBSERVED VARIANCE ■ 9,40695 

PREDICTED MEAN • 2.89981 PREDICTED VARIANCE • 9,27860 


gamma DIST 

PARAMETERS! 

ALPMA ■ 

0.67068 BETA 

■ 0,25864 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVED 

OBSERVED 

expectfd 

DIFFERENCE 

1 

129 

0.4i346 

0,43796 

-0,02450 

2 

65 

0,62179 

0,61528 

0,00651 

3 

40 

0.75000 

0.727D9 

0,02291 

4 

22 

0.82051 

0,80286 

0,01765 

5 

13 

0.86218 

0,85598 

0,00620 

6 

8 

0.88782 

0.89396 

-0,00614 

7 

8 

0,91346 

0.92147 

•0,00801 

R 

4 

0.92628 

0.94159 

-0,01531 

9 

5 

0.94231 

0,95640 

•0,01409 

10 

5 

0.95833 

0,96735 

•0,00902 

11 

2 

0,96474 

0,97550 

-0,01076 

12 

5 

0,98077 

0,98157 

•0,00080 

13 

2 

0.9R718 

0.98612 

0,00106 

14 

1 

0.99038 

0.98952 

0,00086 

15 

1 

0,99359 

0,99209 

0.00150 

16 

0 

0.99359 

0.99401 

-0,00042 

17 

1 

0.99679 

0,99547 

0,00133 

1R 

0 

0.99679 

0,99656 

0,00023 

19 

0 

0.99679 

0,99789 

•0,00060 

20 

0 

0.99679 

0,99802 

•0,00122 

OVER 20 

1 

1,00000 

l.ooooo 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 111 
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KDIMOGOROV-SMIRNOV GOODNESS OF PIT TEST 
USING PREDICTED VAlUFS 


A-LEVEL • 2.500 AUTOCORR E L&T I ON PARAMETER ■ 0.075 


OBSERVED MEAN - 2.12000 OBSERVED VARIANCE ■ 3.30559 

PREDICTED MEAN ■ 2,05026 PREDICTED VARIANCE ■ 2,62239 


GAMMA 0 1 ST 

PARAMETERS* 

ALPHA ■ 

0.91646 BFTA ■ 

0,59116 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

OBSERVED 

OBSERVFO 

pxpectfd 

DIFFERENCE 

1 

25 

0.50000 

0.40842 

0,01150 

2 

14 

0. 78000 

0,72596 

0,05406 

3 

4 

0.86000 

0.85167 

0,00833 

4 

3 

0.92000 

0.91930 

0,00070 

5 

l 

0.94000 

0,95596 

.0, 01596 

6 

1 

0.96000 

0.97591 

*0,01591 

7 

0 

0.96000 

0,90600 

•0,02680 

8 

1 

0.98000 

0,99276 

•0,01276 

9 

0 

0.90000 

0,99663 

-0,01603 

10 

1 

1.00000 

0,99702 

0,00210 

11 

0 

1,00000 

0,99800 

0,00120 

12 

0 

l. ooooo 

0,99934 

0,00066 

13 

0 

1.00000 

0.99964 

0,00036 

U 

0 

1.00000 

0.99900 

0,00620 

15 

0 

i. ooooo 

0,99909 

0,00011 

16 

0 

1.00000 

0,99996 

0,00006 

17 

0 

1.00000 

0,99997 

0,00003 

18 

0 

1.00000 

0,99998 

0,00002 

19 

0 

1.00000 

0,99999 

0,00001 

20 

0 

1.00000 

0,99999 

0,00001 

OVER 20 

0 

1.00000 

1,00000 

0,00000 

TnTAL NUMBER 

OF OVERSHOOTS • 30 
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KOLMOGOAOV-SMIANOV GOODNESS of fit test 
USING PAfOICTEO VALUFS 


A-LEVEL • 0,000 AUTOCORRELATION PARAMETER • 0.250 


OBSERVED MEAN • 4,34698 OBSERVED VARIANCE ■ 20,87042 

PREDICTED MEAN ■ 4,28208 PREDICTED VARIANCE • 21,62711 


gamma DIST 

PARAMETERS! 

ALPHA ■ 

0.66140 BETA 

• 0,17488 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

observed 

observed 

expected 

difference 

1 

890 

0.31415 

0,32685 

-0,01270 

2 

485 

0.48535 

0.48415 

0.00120 

3 

311 

0.59513 

0,59445 

0,00068 

4 

218 

0.6T208 

0.67692 

-0,00484 

5 

187 

0.73809 

0,74044 

-0,00235 

6 

135 

0.78574 

0,79024 

-0,00451 

7 

130 

0.83163 

0,87975 

0,00188 

8 

79 

0.85951 

0,84134 

-0,00182 

9 

74 

0.88563 

0.88675 

-0,00112 

10 

58 

0.90611 

0,90730 

-0,00119 

11 

52 

0.92446 

0,92397 

0,00049 

12 

44 

0.93999 

0,93755 

0,00245 

13 

27 

0.94952 

0,94863 

0,00090 

14 

26 

0.95870 

0,95769 

0,001C1 

15 

26 

0.96788 

0.96511 

0,00277 

16 

17 

0.97388 

0,97121 

0,00267 

17 

9 

0.97706 

0,97622 

0.00084 

18 

13 

0.98164 

0.98034 

0,00131 

19 

7 

0.98412 

0.98373 

0,00038 

20 

10 

0.98765 

0,98654 

0,00111 

OVER 20 

35 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER of OVERSHOOTS ■ 2831 


51 



KnLMOGOROV-SMlRNOV GOODNESS QF FIT TEST 
USING PREDICTED VALUFS 


A-LEVEL ■ 0.750 AUTOCORRELATION PARAMETER • 0.250 

OBSERVED MEAN ■ 2.63016 OBSERVED VARIANCE • 3.7056) 

PREDICTED MEAN - 2.76463 PREDICTED VARIANCE • 7,65495 


GAMMA DIST 

PARAMETERS* 

ALPHA • 

0.66997 BETA ■ 

0,29384 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

OBSERVEO 

OBSERVED 

EXPECTfD 

difference 

1 

929 

0.44049 

0,47646 

0,00403 

2 

423 

0.64106 

0,623*3 

0,01 U<* 

3 

266 

0.76719 

0,74067 

0,02652 

4 

164 

0.04495 

0.01827 

0,02660 

5 

110 

0.89711 

0,87134 

0,02577 

6 

68 

0.92935 

0.90827 

0,02100 

7 

48 

0.95211 

0,97426 

0,01785 

e 

35 

0.96071 

0.93269 

0,01601 

9 

23 

0.97961 

0,96505 

0,01376 

10 

18 

0.90813 

0.97528 

0,01287 

1 1 

7 

0.99147 

0.98207 

0,00940 

12 

6 

0.99431 

0,90697 

0.00734 

1 3 

1 

0.99478 

0.99051 

0,00427 

14 

3 

0.99621 

0.99309 

0,00712 

15 

1 

0.99660 

0,99495 

0,00173 

16 

1 

0.99716 

0.99631 

0,00084 

17 

3 

0.99038 

0.99730 

0,00127 

IB 

0 

0.99058 

0,99803 

0,00055 

19 

1 

0.99905 

0.99855 

0.00050 

20 

0 

0.99903 

0,99894 

0,00011 

OVER 20 

2 

1.00000 

1.00000 

0,00000 

TnTAl NUMBER 

OF OVERSHOOTS • 2109 
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KOLMOGOAOV-SMIRNOV GOODNESS of fit test 
USINO PREDICTED values 


A-LEVEL ■ 1.250 AUTOCDRRELAY ION PARAMETER ■ 0*250 

OBSERVED MEAN ■ 2*0)785 OBSERVED VARIANCE • 2*76502 

PREDICTED MEAN • 2.167*5 PREDICTED VARIANCE • 3,39810 


GAMMA DIST 

PARAMETERS! 

ALPHA ■ 

0.81822 BETA ■ 

0,69070 


NUMBER 

CUMULATIVE 

CUMULATIVE 


OURATIQN 

OBSERVED 

OBSERVED 

expected 

DIFFERENCE 

1 

655 

0,53865 

0.66331 

0,0553* 

2 

25* 

0.7*753 

0.70*69 

0,0*205 

3 

139 

0 . 06 1 8 A 

0.82761 

0,03*2* 

* 

73 

0.92168 

0.69828 

0.02359 

5 

*2 

0.956*1 

0,93958 

0,01683 

6 

22 

0.97*51 

0.96395 

0,01055 

7 

17 

0.960*9 

0.978*2 

0,01007 

8 

3 

0.99095 

0.9R705 

0,00390 

9 

5 

0.99507 

0.99221 

0,00285 

10 

2 

0.99671 

0,99531 

0,00160 

11 

1 

0.99753 

0,99717 

0,00036 

12 

2 

0*99918 

0.99829 

0,00089 

13 

0 

0.99918 

0.99897 

0,00021 

U 

0 

0.99918 

0.99938 

•0,00020 

15 

0 

0*99918 

0.99962 

•0,000*6 

16 

0 

0.99918 

0,99977 

-0,00059 

17 

0 

0.99918 

0,99986 

•0,00068 

18 

0 

0.99918 

0.99992 

-0,0007* 

19 

0 

0.99918 

0.99995 

•0,00077 

20 

0 

0.99918 

0.99997 

-0,00079 

OVER 20 

1 

1.00000 

1.00000 

0,00000 

TOTAL NUMBER 

OF OVERSHOOTS - 1216 
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KnLMnr.OROV-SMlRNOV goodness of fit test 
USING PREDICTED VALUFS 


A-LEVC 

L * 1.750 

AUTOCORRELATION parameter ■ 

0.250 

observed 

mean * 1 

.72134 OBSERVED VARIANCE ■ 

1,43816 

PREDICTED 

MEAN * 1 

.76625 PREDICTED VARIANCE ■ 

1,40658 

gamma DIST 

PARAMETERS! ALPHA ■ 

1.13992 BETA 

• 0,90023 


number 

CUMULATIVE 

CUMULATIVE 


duration 

observed 

observed 

expected 

DIFFERENCE 

1 

307 

0.60672 

0.53131 

0,07541 

2 

115 

0.83399 

0.79747 

0,03653 

3 

41 

0.91502 

0,91418 

0,00004 

4 

22 

0.95850 

0,96400 

*0,00550 

5 

10 

0.97826 

0,98498 

-0,00672 

6 

5 

0.98814 

0.99376 

*0,00562 

7 

5 

0.99802 

0,99742 

0,00061 

8 

1 

1.00000 

0.99893 

0,00107 

9 

0 

1.00000 

0.99956 

0,00044 

10 

0 

1.00000 

0,99982 

0,00018 

U 

0 

1.00000 

0,99993 

0,00007 

12 

0 

1.00000 

0.99997 

0,00003 

13 

0 

1.00000 

0.99999 

0,00001 

14 

0 

1.00000 

0.99999 

0,00001 

15 

0 

1.00000 

r, ooooo 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

l.ooooo 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1.00000 

0,00000 


total number of OVERSHOOTS ■ S06 


I I 


I 
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KfUMOGOROV-SMlRNOV GOODNESS OF FIT TEST 

using predicted values 


A-LEVEL ■ 2.500 AUTOCORRELATION PARAMETER ■ 0.250 

OBSERVED MEAN « 1.38272 OBSERVED VARIANCE ■ 0,85353 

PREDICTED MEAN • 1.39699 PREDICTED VARIANCE • 0,48802 


GAMMA 3 I ST 

PARAMETERS! 

ALPHA ■ 

1,64867 BETA - 

1,83802 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

OBSERVEO 

observed 

expected 

DIFFERENCE 

1 

64 

0.79012 

0.65623 

0.13389 

2 

10 

0.91358 

0.92382 

*0,01024 

3 

3 

0.95062 

0.98494 

-0,03432 

A 

2 

0.97531 

0.99718 

*0,02187 

5 

1 

0.98765 

0.99949 

-0.01184 

6 

1 

1.00000 

0,99991 

0,00009 

7 

0 

1.00000 

0,99998 

0,00002 

8 

0 

1,00000 

1,00000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1.00000 

1,00000 

0, OO'iOO 

11 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

0,00000 

14 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 81 
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KOLMOGOROV- SMIRNOV GOODNESS OF FIT TEST 
USING PREDICTED VAlUFS 


A-LEVEL • 0.000 AUTOCORRELATION PARAMETER ■ 

OBSERVED MEAN • 2.89229 OBSERVED VARIANCE • 

PREDICTED MEAN * 2,77466 PREDICTED VARIANCE • 

GAMMA DIST PARAMETERS* ALPHA ■ 1,04078 BETA 



NUMBER 

CUMULATIVE 

CUMULATIVE 

DURATION 

OBSERVED 

observed 

expectfu 

1 

1628 

0.37228 

0,34775 

2 

912 

0.58084 

0,58128 

3 

610 

0.72033 

0.73232 

4 

412 

0.81454 

0,82925 

5 

280 

0.87857 

0,89123 

6 

161 

0.91539 

0,93078 

7 

122 

0.94329 

0,95598 

8 

93 

0.96456 

0,97202 

9 

43 

0.97439 

0.98222 

10 

39 

0.98331 

0,98871 

11 

20 

0.98788 

0,99283 

12 

14 

0.99108 

0,99545 

13 

12 

0.99383 

0.99711 

14 

10 

0.99611 

0,99817 

15 

8 

0.99794 

0,99884 

16 

2 

0.99840 

0,99926 

17 

2 

0.99886 

0.99953 

18 

1 

0.99909 

0.99970 

19 

0 

0.99909 

0,99981 

20 

1 

0.99931 

0,99988 

OVER 20 

3 

1.00000 

1,00000 

TOTAL NUMBER 

OF OVERSHOOTS ■ 4373 



0,750 

6,06546 

4,97132 

0,45756 


DIFFERENCE 

0,02453 

-0,00045 

-0,01199 

-0,01470 

-0,01265 

-0,01539 

-0,01269 

-0,00746 

-0.00783 

-0.00540 

-0,00495 

-0,00437 

-0,00329 

-0.00206 

-0,00090 

•0,00086 

-0,00068 

-0,00062 

-0,00073 

-0,00057 

0,00000 
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KOLMOGOROV- SMIRNOV GOODNESS OF FIT TEST 
USINO PREDICTED VAlUFS 


A-LEVEL ■ 0.750 AUTOCORRELATION PARAMETER • 0.750 

OBSERVED MEAN « 1.78643 OBSERVED VARIANCE ■ 1,62950 

PREDICTED MEAN • 1,81846 PREDICTED VARIANCE ■ 1,71390 


GAMMA DIST 

PARAMETERS! 

ALPHA ■ 

1.01425 BETA 

■ 0,76927 


number 

CUMULATIVE 

CUMULATIVE 


duration 

observeo 

OBSERVED 

expected 

DIFFERENCE 

1 

1860 

0.58161 

0.52991 

0,05171 

2 

743 

0.81395 

0.78086 

0,03308 

3 

305 

0.90932 

0.89804 

0,01127 

4 

141 

0.95341 

0.95261 

0,00080 

5 

75 

0.97686 

0,97799 

-0,00113 

6 

39 

0.98906 

0.98978 

-0,00072 

7 

20 

0.99531 

0,99525 

0,00006 

8 

7 

0.99750 

0,99780 

•0,00030 

9 

3 

0.99844 

0.99898 

-0,00054 

10 

2 

0.99906 

0.99953 

-0,00046 

11 

1 

0.99937 

0.99978 

-0,00041 

12 

1 

0.99969 

0.99990 

-0,00021 

13 

1 

1.00000 

0,99995 

0,00005 

14 

0 

1.00000 

0,99998 

0.00002 

15 

0 

1.00000 

0,99999 

0,00001 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 3198 
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KOLMOGOROV-SMIRNOV GOODNESS QF FIT TEST 
USING PREDICTED VALUES 


A-LFVEL * 1.250 AUTOCORRELATION PARAMETER ■ 

OBSERVED MEAN * 1.45505 OBSERVED VARIANCE ■ 

PREDICTED MEAN « 1.47133 PREDICTED VARIANCE ■ 

GAMMA DIST PARAMETERS! ALPHA ■ 1,28920 BETA ■ 


0.750 

0.68482 

0,73103 

1,32726 


duration 

NUMBER 

OBSERVEO 

CUMULATIVE 

observed 

CUMULATIVE 

EXPECTED 

difference 

1 

2 

1254 

0.69167 

0,62912 

0,06256 

383 

0.90292 

0 . 8*694 

0,01598 

3 

116 

0.96691 

0,96711 

- 0,00021 

4 

38 

0.9B787 

0,99065 

-0,00278 

5 

16 

0.99669 

0,99738 

-0,00069 

A 

4 

0.99B90 

0.99927 

-0,00037 

u 

7 

2 

1.00000 

0.999B0 

0,00020 

g 

0 

1.00000 

0,99994 

0,00006 

9 

0 

1.00000 

0,99990 

0,00002 

10 

0 

1,00000 

1,00000 

0,00000 

1 1 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

1 3 

0 

1.00000 

1,00000 

0,00000 

U 

0 

1.00000 

1.00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

l, ooooo 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 

TOTAL NUMBER 

OF OVERSHOOTS ■ IBIS 




58 



KHIMOGOROV-SMIRNOV GOOONeSS OF FIT TEST 
USINO predicted valufs 


A-LEVEL » 1.750 AUTOCORRELATION PARAMETER ■ 0.750 

OBSERVED mean * 1 *29041 OBSERVED VARIANCE ■ 0,39792 

PREDICTED MEAN « 1.25889 PREDICTED VARIANCE ■ 0,28063 


gamma D!ST 

parameters 1 

ALPHA ■ 

2.05222 BETA • 

2,70423 


number 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVED 

OBSERVED 

EXPECTfD 

difference 

1 

595 

0.78187 

0.74012 

0,04175 

2 

125 

0.94612 

0.94907 

- 0,02295 

3 

30 

0.98555 

0,99700 

- 0 , C 1146 

4 

9 

0.99737 

0.99974 

- 0,00236 

5 

1 

0.99869 

0,99998 

- 0,00129 

6 

1 

1.00000 

1,00000 

0,00000 

7 

0 

1.00000 

1,00000 

0,00000 

8 

0 

1.00000 

l ,00000 

O.OOOQO 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

11 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

l'.ooooo 

0,00000 

14 

0 

1.00000 

l.ooooo 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

l . 00000 

1,00000 

0,00000 

19 

0 

1,00000 

1.00000 

0 , C' ooo 

20 

0 

1.00000 

1,00000 

o.ocooo 

OVER 20 

0 

1,00000 

1,00000 

0,00000 


TOTH NUMBER OF OVERSHOOTS ■ 761 
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KHLMOGORDV-SMlRNOV GOODNESS OF FIT TEST 
USING PREDICTED VALUES 


A-LEVEL « 2.500 AUTOCORRELATION PARAMETER ■ 0.750 

OBSERVED MEAN « 1.10*89 OBSERVED VARIANCE ■ 0,1*984 

PREDICTED MEAN - 1 • 10*43 PREDICTED VARIANCE ■ 0,07902 


GAMMA D I ST 

PARAMETERS! 

ALPHA ■ 

4,65376 BETA ■ 

7,67*02 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

observed 

observed 

EXPECTED 

DIFFERENCE 

l 

131 

0.91608 

0,90746 

0.00A62 

2 

10 

0.9B601 

0.99959 

-0,01357 

3 

1 

0.99301 

1,00000 

-0,00699 

* 

1 

1.00000 

1.00000 

0,00000 

5 

0 

1.00000 

1,00000 

0,00000 

6 

0 

1.00000 

1,000^0 

0,00000 

7 

0 

1.00000 

1,00000 

0,00000 

fl 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1.00000 

l.ooooo 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

11 

0 

1.00000 

l.ooooo 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

0,00000 

U 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

l.ooooo 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1,00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 149 
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KOIMOOOROV-SMIRNOV GOODNESS OF FIT TEST 
USINO PREDICTED VAlUFS 


A-LEVEL - 0.000 AUTOCORRELATION PAR AMFTER ■ 1,000 

OBSERVED MEAN • 2.60607 OBSERVED VARIANCE ■ 4,48415 

PREDICTED MEAN - 2. 5**1* PREDICTED VARIANCE ■ 3,55638 


GAMMA 01 ST 

PARAMETERS! 

ALPHA • 

1.15199 BETA - 

0,56912 


number 

cumulative 

CUMULATIVE 


DURATION 

observed 

observed 

fxpectfd 

difference 

1 

1925 

0.39765 

0.36326 

0,03438 

2 

1104 

0.62370 

0,61857 

0,00713 

3 

689 

0.76802 

0.77530 

-0,00728 

4 

412 

0.85313 

0.86879 

-0,01567 

5 

262 

0.90725 

0.92381 

-0,01656 

6 

178 

0.94402 

0.95592 

-0,01190 

7 

101 

0.96488 

0,97457 

-0,00969 

B 

59 

0.97707 

0.98536 

-0,00829 

9 

43 

0.98595 

0,99159 

-0,00564 

10 

18 

0.98967 

0,99517 

-0,00550 

11 

18 

0.99339 

0.99723 

-0,00384 

12 

9 

0.99525 

0,99842 

-0,00317 

13 

5 

0.99628 

0.99909 

-0,00281 

14 

7 

0.99773 

0,99948 

-0,00175 

15 

5 

0.99876 

0,99970 

•0,00094 

16 

3 

0.99938 

0.99983 

-0,00045 

17 

2 

0.99979 

0.99990 

•0,00011 

IB 

1 

1.00000 

0.99994 

0,00006 

19 

0 

1.00000 

0.99997 

0,00003 

20 

0 

1 .00000 

0,99998 

0,00002 

over 20 

0 

1,00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 4961 
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KflLMOGOROV-SMIRNOV GOODNESS OF FIT TEST 

USING predicted valufs 


A-LEVEL • 0.750 AUTOCOR 0 E l AT I ON PARAMETER ■ 1,000 

OBSERVED MEAN * 1,66955 OBSERVED VARIANCE ■ 1,30336 

PREDICTED MEAN « 1,66799 PREDICTED VARIANCE ■ 1,25*77 


GAMMA 01 ST 

PARAMETERS! 

ALPHA ■ 

1,08720 BETA ■ 

0 , 9308 * 


NUMBER 

CUMULATIVE 

CUMULATIVE 


OURATION 

observed 

observed 

EXPECTED 

DIFFERENCE 

1 

2151 

0.61971 

0.56719 

0,05252 

2 

770 

0 . 8 * 15 * 

0,82263 

0,01891 

3 

302 

0.92855 

0.92822 

0,00033 

* 

1*0 

0.96889 

0 . 9711 * 

- 0,00225 

5 

61 

0 . 986*6 

0 , 988 ** 

- 0,00198 

6 

21 

0.99251 

0.99538 

- 0,00287 

7 

15 

0.99683 

0,99816 

- 0,00133 

8 

7 

0.99885 

0,99927 

- 0 , 000*2 

9 

1 

0 . 9991 * 

0,99971 

- 0,00057 

10 

1 

0 . 999*2 

0.99988 

- 0 , 000*6 

11 

0 

0 . 999*2 

0,99995 

- 0,00053 

12 

0 

0 . 999*2 

0.99998 

- 0,00056 

13 

1 

0.99971 

0 99999 

- 0,00028 

1 * 

1 

1.00000 

1, DODO 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

r . ooooo 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 

total number 

OF OVERSHOOTS ■ 3*71 
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KniMOGOROV-SMlRNOV GOODNESS OF FIT TEST 
USING PREDICTED VALUES 


A-LEVEl « 1,230 AUTOCORRELATION PARAMETER • 1,000 


OBSERVED MEAN » 1,37456 OBSERVED VARIANCE ■ 0,57043 

PREDICTED MIAN ■ 1,36199 PREDICTED VARIANCE ■ 0,54155 


GAMMA 01 ST 

parameters l 

ALPHA ■ 

1 , 37265 BFTA 

■ 1,59172 


number 

CUMULATIVE 

CJMI'LATJVE 


duration 

OBSERVED 

observed 

fxpectfd 

DIFFERENCE 

1 

1485 

0.74065 

0.67827 

0,06238 

2 

363 

0.92170 

0.92107 

0,00063 

3 

105 

0,97406 

0,98184 

-0,00777 

4 

37 

0.99252 

0,99595 

-0,00343 

5 

10 

0.99751 

0,99911 

-0,00161 

6 

4 

0.99950 

0,99981 

-0,00031 

7 

0 

0.99950 

0,99996 

-0 , 00046 

8 

1 

1.00000 

0,99999 

0,00001 

9 

0 

1.00000 

i.ooooo 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

U 

0 

1,00000 

1,00000 

0,00000 

12 

0 

1.00000 

i.ooooo 

0,00000 

13 

0 

1,00000 

1.00000 

0,00000 

14 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

i.ooooo 

1,00000 

0,00000 

17 

0 

1,00000 

I.OOOOO 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1,00000 

1,00000 

0,00000 

20 

0 

1,00000 

1,00000 

0,00000 

OVER 20 

0 

1,00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS - 2005 
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KIUMOGOROV-SMIRNOV GQOONESS OF FIT TEST 
USING predicted values 


A-LEVEL ■ 1.750 AUTOCORRELATION PARAMETER • 1,000 

OBSERVED MEAN « 1,22261 OBSERVED VARIANCE ■ 0,32530 

PREDICTED MEAN ■ 1,17950 PREDICTED VARIANCE • 0,20484 


GAMMA DIST 

PARAMETERS! 

ALPHA • 

2.25402 BETA 

• 3,31716 


number 

CUMULATIVE 

CUMULATIVE 


duration 

observeo 

observed 

expectfd 

DIFFERENCE 

1 

722 

0.83276 

0.79859 

0,03417 

2 

no 

0.95963 

0,95515 

.•0,02552 

3 

26 

0.98962 

0.99915 

-0,00954 

4 

6 

0.99654 

0,99996 

-0,00342 

5 

2 

0,99085 

i.onooo 

-0,00115 

6 

1 

1.00000 

l.ooooo 

0,00000 

7 

0 

1.00000 

1,00000 

0,00000 

8 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1,00000 

1,00000 

0,00000 

11 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1,00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

0,00000 

u 

0 

1,00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1,00000 

1.00000 

0,00000 

TOTAL NUMBER 

OF OVERSHOOTS ■ 867 
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KniMOGOROV-SHlRNUV GOODNESS 0* FIT TEST 
US1N0 PREDICTED VALUES 


A-LEVEL ■ 2.500 AUTOCORRELATION PARAMETER • 1,000 


OBSERVED mean • 1 . 1 3043 OBSERVE* VARIANCE • 0.1A2A1 

PREDICTED MEAN « 1 .061 AO PREDICTED VARIANCE ■ 0.0A707 


GAMMA DtST 

PARAMETERS! 

ALPHA ■ 

6.69515 BETA 

• 11,92585 


number 

CUMULATIVE 

CUMULATIVE 


OUR AT ION 

observed 

OBSERVED 

expected 

DIFFERENCE 

i 

122 

0.8BA06 

0.96205 

-0,07799 

2 

1A 

0.98551 

0.99999 

•0.01AA8 

3 

2 

1.00000 

1,00000 

0,00000 

A 

0 

l. 00000 

1,00000 

0,00000 

5 

0 

1.00000 

1,00000 

0,00000 

6 

0 

1.00000 

1,00000 

0,00000 

7 

0 

1.00000 

1,00000 

0,00000 

8 

0 

1.00000 

1,0*000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

U 

0 

l.ooooo 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

o.ooooo 

1A 

0 

1.00000 

l.ooooo 

o.ooooo 

13 

0 

1.00000 

1, 00070 

0,00000 

16 

0 

1.00000 

l.OOO^o 

0,00000 

17 

0 

1.00000 

1.00000 

0,00000 

IB 

0 

1.00000 

1,00000 

0,00000 

19 

0 

l.ooooo 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 1)1 
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KntMOGOROV-SMlRNQV GOOONeSS OF FIT TEST 

USING predicted valufs 


A-LEVEL • 0.000 AUTOCORRELATION PARAMETER ■ 2.500 


OBSERVED MEAN - 2.07269 OBSERVED VARIANCE * 2,21012 

PREDICTED MEAN - 2.01*08 PREDICTED VARIANCE • 1,82638 


gamma oist 

PARAMETERS! 

ALPHA ■ 

1,26181 BETA • 

0,83119 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVED 

OBSERVED 

expected 

DIFFERENCE 

1 

2890 

0. *8295 

0 . **830 

0,03*65 

2 

1526 

0.73797 

0.73182 

0,0061* 

3 

703 

0.855*5 

0.87*03 

-0,01858 

* 

*19 

0.925*7 

0.9*186 

-0,01639 

5 

226 

0.9632* 

0.973*5 

-0,01022 

6 

102 

0.98028 

0.9*797 

•0,00769 

7 

60 

0.99031 

0,99*57 

•0,00*27 

8 

3* 

0.99599 

0.99756 

-0,00157 

9 

9 

0.997*9 

0.99891 

-0,001*2 

10 

8 

0.99883 

0.99951 

•0,00068 

11 

* 

0.99950 

0.99978 

•0,00028 

12 

2 

0.99983 

0.99990 

-0,00007 

13 

1 

1.00000 

0.99996 

0.O0O0* 

1* 

0 

1.00000 

0.99998 

0,00002 

15 

0 

1.00000 

0,99999 

0,00001 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1.00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

l. 00000 

0,00000 

20 

0 

1.00000 

1 ,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 598* 



KHtMOGOBOV-SMUNOV GUODNESS OF FIT TEST 

USING predicted values 


A-LEVEl • 0*750 AUTOCORRELATION PARAMETER • 2,300 

OBSERVED mean • 1,33962 OBSERVED VARIANCE ■ 0,44402 

PREDICTED MEAN • 1,30232 PREDICTED VARIANCE • 0, 6*758 


gamma dIST 

PARAMETERS! 

ALPHA • 

1.20215 BETA * 

1,36249 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVED 

OBSERVED 

PXPECTFD 

DIFFERENT.! 

1 

3204 

0,7*581 

0.67181 

0,07400 

2 

913 

0,93506 

0.90726 

0,02780 

3 

206 

0.98301 

0,97446 

0,00835 

* 

60 

0,99697 

0.99319 

0,00379 

5 

11 

0.99953 

0.99019 

0,00135 

6 

2 

1.00000 

0.99952 

0,000^6 

7 

0 

1.00000 

0.99987 

0,00013 

8 

0 

l .00000 

0,99997 

0,00003 

9 

0 

1,00000 

0.99999 

0,00001 

10 

0 

1.00000 

1,00000 

0,00000 

1 1 

0 

1.00000 

l.OOOOp 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1.00000 

0,00000 

1* 

0 

1.00000 

1,00000 

0 • 00000 

13 

0 

1.00000 

1.00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1,00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1.00000 

0,00000 

20 

0 

1,00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

l.ooooo 

0,00000 


TOTAL NUMBER OP 0VERJM00T1 • 4296 
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KOLMOGOROV-SMlRNOV GOODNESS OR PIT TEST 
USING PREDICTS VALUES 


A-LEVEL * 1.250 AUTOCORRELATION PARAMETER » 2,500 

OBSERVED MEAN ■ 1 , 14412 OBSERVED VARIANCE ■ 0,15547 

PREDICTED MEAN • 1,10076 PREDICTED VARIANCE ■ 0,27002 


GAMMA D! ST 

PARAMETERS! 

ALPHA > 

1,61690 BETA 

■ 2,44703 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVED 

OBSERVED 

EXPECT FD 

DIFFERENCE 

1 

2061 

0.B7109 

0,79418 

0,07671 

2 

271 

0.98563 

0,97569 

0,01054 

3 

32 

0.99915 

0,99732 

0,00183 

4 

2 

1.00000 

0,99973 

0,00027 

5 

0 

1.00000 

0,99997 

0,00003 

6 

0 

1,00000 

1.00000 

0,00000 

7 

0 

1.00000 

1,00000 

0,00000 

8 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1.00000 

l.ooooo 

0,00000 

11 

0 

1,00000 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

o.ooooo 

13 

0 

1.00000 

1,00000 

0,00000 

14 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

l.ooooo 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1,00000 

1,00000 

0,00000 

OVER 20 

0 

1,00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 2366 
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KOLMOGOROV-SMIRNOV GOODNESS of fit test 
USING PREDICTED VALUrS 


A-LEVEL « 1.730 AUTOCORRELATION PARAMETER ■ 2.500 

OBSERVED MEAN - 1.03657 OBSERVED VARIANCE ■ 0,05339 

PREDICTED MEAN - 1.03711 PREDICTED VARIANCE • 0,08885 


gamma oist 

PARAMETERS) 

ALPHA • 

3,24706 BETA ■ 

6.04538 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

OBSERVED 

OBSERVFD 

expectfd 

DIFFERENCE 

1 

935 

0.94444 

0.92290 

0,02134 

2 

54 

0.99899 

0.99928 

-0,00029 

3 

1 

1.00000 

1 . OnOfiQ 

0,00000 

4 

0 

1.00000 

UOOOOO 

o , ooooo 

5 

0 

1.00000 

l.ODODO 

0,00000 

6 

0 

1.00000 

1,00000 

0,00000 

7 

0 

1.00000 

1,00000 

0,00000 

8 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1,00000 

1,00000 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

11 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1.00000 

l.ooooo 

0,00000 

13 

0 

1.00000 

1,00000 

0,00000 

14 

0 

1.00000 

1.00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

l.ooooo 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 990 
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KOLHOGOROV-SHIRNQV GOODNESS of fit test 
USING PREDICTED VALUES 


A-LEVEl ■ 2.500 AUTOCORRELATION PARAMETER ■ 2,500 

OBSERVED MEAN ■ 1,01527 OBSERVED VARIANCE ■ 0,0150* 

PREDICTED MEAN ■ 0.98398 PREDICTED VARIANCE • 0,00719 


gamma 0 1 ST 

PARAMETERS! 

ALPHA ■ 

32.56425 BETA ■ 

67,28482 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVED 

OBSERVED 

expected 

OIFFERENCl 

1 

129 

0.98473 

1.00000 

-0,01527 

2 

2 

1.00000 

1,00000 

0,00000 

3 

0 

1.00000 

1.00000 

0.00000 

A 

0 

1.00000 

1,00000 

0,00000 

5 

0 

1.00000 

1,00000 

0 , ooooo 

6 

0 

1.00000 

1,00000 

0,00000 

7 

0 

1.00000 

1,00000 

0,00000 

8 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

l.ooooo 

1,00000 

0,00000 

11 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

0.00000 

14 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00001 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 1J1 
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KOLMOGOROV-SMIRNQV GOODNESS QF FIT TEST 
USING PREDICTED VALUFS 


A-LEVEL ■ 0.000 AUTOCORRELATION PARAMETER • 3,000 


OBSERVED MEAN • 1.94846 OBSERVED VARIANCE ■ 1,76810 

PREDICTED MEAN ■ 1,88194 PREDICTED VARIANCE ■ 1,72155 


GAMMA DIST 

PARAMETERS! 

ALPHA » 

1,10933 BETA ■ 

0,80273 


NUMBER 

cumulative 

CUMULATIVE 


duration 

OBSERVED 

OBSERVED 

EXPECTFD 

DIFFERENCE 

1 

3218 

0.50877 

0.50156 

0,00722 

2 

1568 

0.75984 

0,76606 

-0.00622 

3 

791 

0.88490 

0,89180 

-0,00690 

4 

384 

0.94561 

0.95033 

-0,00471 

5 

180 

0.97407 

0,97730 

-0,00323 

6 

87 

0.98783 

0.98966 

-0,00183 

7 

37 

0.99368 

0.99530 

-0,00162 

8 

26 

0.99779 

0,99787 

-o.ooooe 

9 

10 

0.99937 

0,99903 

0,00033 

10 

3 

0.99984 

0.99956 

0,00028 

11 

0 

0.999B4 

0,99980 

0,00004 

12 

1 

1.00000 

0.99991 

0,00009 

13 

0 

1.00000 

0,99996 

0,00004 

14 

0 

1.00000 

0.99998 

0,00002 

15 

0 

1.00000 

0.99994 

0,00001 

16 

0 

1.00000 

1,OOOOC 

0,00000 

17 

0 

1,00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1.00000 

0,00000 

20 

0 

1.00000 

1.00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 6325 
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KDIMOGOROV-SMIRNQV GOODNESS OP PIT test 
USING PREDICTED VALUES 


A-LEVEL = 0.750 AUTOCORRELATION PARAMETER ■ 5,000 

OBSERVED MEAN * 1. 294*6 OBSERVED VARIANCE ■ 0,37522 

PREDICTED MEAN ■ 1.33*85 PREDICTED VARIANCE ■ 0,42310 


gamma oist 

parameters 1 

ALPHA ■ 

1,6*730 BETA * 

1.97316 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

observed 

observed 

expected 

difference 

1 

3348 

0.77019 

0.69057 

0,07961 

2 

786 

0.95100 

0.93979 

0,01121 

3 

159 

0.98758 

0,9*958 

- 0,00200 

4 

*1 

0.99701 

0,99829 

* 0,00128 

5 

12 

0.99977 

0.99973 

0,00004 

6 

1 

1.00000 

0.99996 

0,00004 

7 

0 

1.00000 

0.99999 

0,00001 

B 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1,00000 

1.00000 

0,00000 

11 

0 

1.00000 

1.00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

0,00000 

1* 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ *3*7 
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KOLMOGOROV- SMIRNOV GOODNESS OF FIT TEST 
USING PREDATED VALUES 


A-LEVEL ■ 1.250 AUTOCORRELATION PARAMETER • 5,000 

OBSERVED MEAN ■ 1,12570 OBSERVEO VARIANCE ■ 0,13993 

PREDICTED MEAN ■ 1,14004 PREDICTED VARIANCE • 0,11490 


GAMMA DIST 

PARAMETERS! 

ALPHA ■ 

2,56518 beta 

• 5,57021 


NUMBER 

CUMULATIVE 

CUMULATIVE 


duration 

observed 

observed 

expected 

difference 

1 

2066 

0.88631 

0.8401 1 

0,02621 

2 

242 

0.99013 

0.99751 

-0,00738 

3 

20 

0.99871 

0,99998 

-0 ,00126 

4 

1 

0.99914 

l.oocoo 

-0,00086 

5 

2 

1.00000 

1,00000 

0,00000 

6 

0 

1,00000 

1,00000 

0,00000 

7 

0 

1,00000 

1,00000 

0,00000 

B 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1,00000 

1,00000 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

11 

0 

1,00000 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

0,00000 

14 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

1,00000 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1,00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS • 2)21 
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KHIMOGOROV-SHIRNOV GQHDNeSS of FIT test 
using predicted values 


A-LEVEL - 1.750 AUTOCORRELATION PARAMETER ■ 5,000 

nBSERVEO MEAN - 1, 03747 OBSERVED VARIANCE • 0,03821 

PREDICTED MEAN - 1,03076 PREDICTED VARIANCE • 0,01720 


gamma dist 

PARAMETERS! 

ALPHA • 

16.37617 BETA • 

30,85400 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

observed 

observed 

EXPECTED 

DIFFERENCE 

1 

900 

0.96360 

0.99837 

- 0,03478 

2 

33 

0.99893 

1,00000 

* 0,00107 

3 

1 

l.ooooo 

1,00000 

0,00000 

4 

0 

1.00000 

1,00000 

0,00000 

5 

0 

1.00000 

1,00000 

0,00000 

6 

0 

1.00000 

1.00000 

0,00000 

7 

0 

1.00000 

1,00000 

o.ooooo 

B 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

11 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1.00000 

l.ooooo 

0,00000 

13 

0 

1.00000 

l.ooooo 

0,00000 

14 

0 

1.00000 

1,00000 

0,00000 

15 

0 

1.00000 

l.ooooo 

0,00000 

16 

0 

1.00000 

1,00000 

0,00000 

17 

0 

1,00000 

1,00000 

0,00000 

18 

0 

1.00000 

1,00000 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1.00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOTS ■ 934 
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KfllMOGOROV-SMIRNOV GODDESS QF FIT TEST 
USING PREDICTED VALUES 


A-LEVEL » 2.500 AUTOCORRELATION PARAMETER • 5,000 

OBSERVED MEAN ■ 1,0070* OBSERVEO VARIANCE ■ 0.00699 

PREDICTED MEAN > 0.98610 PREOICTEO VARIANCE ■ 0,00019 


GAMMA D t ST 

PARAMETERS! 

ALPHA *1238 

,77160 BETA 

■2568 , 396<»0 


NUMBER 

CUMULATIVE 

CUMULATIVE 


DURATION 

OBSERVEO 

observed 

EXPECTED 

DIFFERENCE 

1 

161 

0.99296 

l.ooooo 

-0.00706 

2 

1 

l. 00000 

1,00000 

0,00000 

3 

0 

1.00000 

1,00000 

0,00000 

* 

0 

1.00000 

1,00000 

0,00000 

5 

0 

1.00000 

1, 0O0O0 

0,00000 

6 

0 

l. ooooo 

l.ooooo 

0 , OOOOO 

7 

0 

1.00000 

1,00000 

0,00000 

6 

0 

1.00000 

1,00000 

0,00000 

9 

0 

1.00000 

1,00000 

0,00000 

10 

0 

1.00000 

1,00000 

0,00000 

11 

0 

1.00000 

1,00000 

0,00000 

12 

0 

1.00000 

1,00000 

0,00000 

13 

0 

1.00000 

1,00000 

0,00000 

U 

0 

1.00000 

l.ooooo 

0,00000 

15 

0 

l.ooooo 

l, ooooo 

0,00000 

16 

0 

1.00000 

1.00000 

0,00000 

17 

0 

1.00000 

1,00000 

0,00000 

ie 

0 

1.00000 

r, ooooo 

0,00000 

19 

0 

1.00000 

1,00000 

0,00000 

20 

0 

1.00000 

1,00000 

0,00000 

OVER 20 

0 

1,00000 

1,00000 

0,00000 


TOTAL NUMBER OF OVERSHOOT! • 161 
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APPENDIX C 
PREDICTION PROGRAM 
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NOTE DO NOT USE THIS FORM 
WHEN ORDERING NEW OR REVISED 
CARO FORM COMPOSITION 
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«n E 

o» ?; 
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«n a 

cn a 
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cn £ 

cn p 

on 7 

cn E 

«n S 

en r: 

on p 
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on 

cn £ 
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cn P 

on P 

cn f 

cn P 

on f 

on T 

cn a 

(D a 
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on a 
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on z 

on Z 
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cn C 

cn P 
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on a 
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cn r 

cn r 

on r 

on p 
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on p 
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on r 

on = 
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en ® 

on p 
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on •» 
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On •«* 
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on - 
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«n 
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on 
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r 

( 

C. 

l 

L 

l 

C 

L 

l 


L 

c 

L 

C 


i.KAfi pnprnsE i for a r.ivrN process, caic^latf tht PRnB ah i l i t v of 

l*CFF.nilir, A S { l. F. C T t P CPrsSlMG LfVFl Fnk X TIME UNITS 

HI- yUi I- f |.|Pl'T I I) f« iu bS M fc A if MU 

i) PkfiCFSS ^ r A 4 H A j () DlVTATinri, S Tr*£ V 
3 ) AiJTJCnwRFLAn.JM PAPAM^TtO/ BFTA 
A) CROSSING LfcVFl, l 
5) DURATION# X 

K T A I MU, I, 

I M TFGP k C»p 

I /Li ASS ir/'Mf uTS 
C - 5 
P « * 

I'Ti'T P M> Af-r TCRS 

A * ') Mil, sT'. r ;,BF IA,l,X 

S IG ’A « SMf V**,? 


L 

C 


l. 

L 


C 

c 


l 

L 

c 

c 

c 

c 


l 


C A L " L A T p A'\)"5Tt n CPMSSIMG LEVEL 

a » v i ■ n -mu) /stl*fv 

Al Vi ■ A.'S ( Al. Ve ) 


CALC .11 A I I C " PF If. IFMTS 
H ■ Al 1 1 • ( * t T A ) 

AO « <..'i? s >no _ 0.31 1 29*R 

Al « , a A 9 AO ♦ O.OPBF l*R 

A 7 • u. 117*3 ♦ 0,017 A 9*R 
,V> , i .r'S ni 7 - 1,0?3*A*« 
Ml « -1 . 1 * 0. 

■3? ■ 0.? ipA - 0. 16723*R 


♦ 

♦ 

r 

•» 

♦ 


O.OAkl 7*k*«? 
0.0?00?*R**2 
U.0101 

J.5A /96*k+*2 
U.3R t( 2*«**2 
U.076A9*k**2 


♦ 0, 00^8 R*R **i 

- 0,n002**R* v 3 

- O,n0131*P**j 

- 0.09S3S*p*-u3 
» 0 , P 1 AqP *9 * »■ 3 

* 0,01A33»B**3 


CALCULATE P°E n K-TFP Mr AN AND VARIANCE 
k « Ati ♦ A 1 * A L V L ♦ A2*ALVL**2 
p-tan ■ e x p ( r ) 

B ■ 'i n * 0 l*Al VL ♦ R2*ALVL**2 

P » ' A k ■ t-t* * 2 


CAtC 'iATf ('.■ ipability or exchtuiuO l pup x m m e units 

Pf-M'i a 1, - f.AMX ( y f PMFAN, P VAR ) 

•J ‘ 1 T P IT • 1 1 ; ' T I ‘ ; C 

••'MTF <P,9S ) BETA,L,ALVL,MU,P IF AN, SIGMA, P VAR, X, PROP 


hMR-ats 


FI’h i AT ( SF 1 0 , 3 ) 

95 F ik *AT( • 1 »,AOX, 'AUTOCORPtUTin .) PAHAMfTfR ■ ' , F 1 1 , 3//A J X, 'CROSSING 
l L f V r L ■ 1 j F 1 1 , 3, • AOJUSTFU LF VF L • » , F 1 1 , 3/ / A 1 X, ' PROG E SS 

PMF A*. ■'»F11,3iI DISTRIBUTION MEAN ■'/F11,3//S1X, •PROCESS 

3 w A n I Ah f F • ' J F 1 1 • 3# ' DISTRIBUTION VARIANCF ■ » » F 1 1 , 3/ / /SOX# • PRPBA 

ABILITY OF EXCEEDING TH£ CROSSING LEVEL FOR • /5AX, F 1 1 , 3, ' TIME UNITS 
. F A . 3 ) 
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C 'ir'Ki.i I VI ci Sir-rj 1 NCI1MPLE1 E GA U ^ 1 A, GAMY ( X # QP ) # — PIJNCHTO Of* 026 
> ' if pmcisioh njHCTinN ga,'(X(T,xbar,vak) 

1 n ! r I T Kf AL*M D ) 

p 'U U PM C ISlQfi VURMX 

nf s v u AC f *? , / VAR 

■'* • ( 1 ■*X f ' AP ) /VAR 

PY c ! >x 

,)t » *'F 

! ' S 1 1 * ' ■ O.PQ 

ir ( >X.C.T,o,no) GO Tn 2 
(Mu/ s '.no 
(.n T | yo 

? IP ( f'F .r.T.'VH)) GO TO A 
ijAf'X * l.no 

G 1 r.J vn 

H ( (pr ,r.f .2O0.D0) .AOO. (DV.LE.HF) ) CD TQ 21 

ir ( h)F. 0 i. 2.00) .AMO. (PY.Gf, DF^3'.00*rS0RT(0F) )) Gn TU AO 

,:*I = nr 

T'r 3 'Ml*nl IIG(OY) - PY - 'Jl&GM ( OA I ♦ 1,00) 

16 ir(nor.Lb.-Po.uo) on to 10 
pro ■ prxP(nuF) 

1 2 r ’f if » I'P G 

n<ii*« . n S ii'; 4 prr, 

' r ,, 3 pr !,*|W / ( pA I ♦ l.DO) 
n'* I » 'I A l ♦ 1,1)0 
! i- ('AT , ■*! .200,00 ) GO T[1 ?5 
I r ( ,v r G , I. T , :)I h) GO TO 13 
|F C'M .01 .?00. DO) GO T [) 2 5 
Ml TO IP 
IT nrn ■ iff. 

If (n h ( /os<JM,LE,l,n-lA) GO TO 15 
t n i; 

10 ’»A| , IMI ♦ l.DO 

ir I'M ,-.I .2^0.00) GO To 25 
Mr « '•’m ♦ in nr.coY/DAi » 

Tj l* 

21 .111* ’),'>( V)f 

70 P» . ■ ( (•Y/MWfcC. 333333333333333 - l.DO ♦ 1 , PO/DH ) *D5QRT ( DH ) 
•1‘I-ITX * Y ; 1 R i ’ X ( D T N ) 

:.Ou ■ n-piix ♦ DSOM 
r,i’ T i ) on 

15 r.v* ■ ''s'im 

' T,J ') ■) 

?S ii'i » o , ,'o*l)AI 
nr • PA I 

' P- T I | p*> 

40 ,)> x « I r * n LPG( DY ) - DY - PlGGM(PF) 

L » f' . J o 
» 0 • 1 » 1 .'VO 

l ■ l.t'O 
D r *n ■ PY 
•)(*»■ ■l.DO 
J 0 K n • 0 Y 

ok ■ i ,no 
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i A k 

n 

nv - »>r 


1 1 a k r* 

r 

I’K 


n 4 l 

a 

m tn AH 

♦ oak*oal 

l"'l 

a 

i* Ik *f)RH 

♦ DAK *()Rt 

I'AM 

« 

I'u^fOAl 

♦ oakp*dah 

!)'«'( 

a 

ir»KP*i'RL 

♦ OAK P*DRH 

I ,r A 

■ 

o A L / 0 R L 


l)f \\ 

a 

dah/orm 



jr (nrr .ig.o.Do) on To 49 
IF('>AnSUnFA-rFH)/DPB) t lEtliD-i3> 00 TO M 

[>K ■ "K ♦ 1 , DO 

r,n rn <*? 

4 i or x ■ r»rx ♦ DloG(OFR) 

GAMX > 1.00 

I F ( OFX .r,£ ,-flO.DO) GAHX • 1,00 m DEXPtQPX) 
on T[1 99 

45 GAMX ■ 1,00 

r.n rn 99 

99 RfTURM 
EMD 
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t.nitir I’MCISIIJN FUNCTION VnRHX(UZ) 

I ' ' n 1 f I T u £ A L * B t D > 

PP t a , V)P9^22RO401 / .?'' 

!) < s t’A' 1 S ( DZ ) 

ir (pv t r. f , *,no> r.n rn in 

n A I = p.l-n 

nt’ i = i . ;to 

;|AM a X 

;i"m = l.r.n 

’.lA'i = w.l'" 

s L)A-J e pA;i ♦ 1 • [>r* 

) a i < -(?.nn*r>AN - l.OQ) *t)X*DX 

nr i s 4,i»n*nAM - 1,00 

1) A L = i:*M*nAH + I) A I * 0 A L 

n n L * r>nj*()RH ♦ l)AI*0PL 

r>A1 = - [)A I 

O'M = 2.0 n ♦ OP! 

■)M s ;>Oi*nAL ♦ D A I * l) A H 
ptM = f*n 1 ♦f)BL + D A I * *.) P H 

liFA = "AL/DPL 

riFR = 1 1 Am / !)i'H 

I r ( m " . r j . u • DP ) &o TU 20 

lF(''.\ uc >((nFR-OFA)/0FA),LE.l.U“l4) GO TQ 2 n 
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1 0 

'ML 

r 

m . i n 


D r ' L 

= 

l’.r o 


CAm 

2 

1 .oo 


orm 

x 

iiX 


D 1 ' I 

z 

nv 


PA , 4 

2 

1 .DO 


nr a 

2 

1 .no. x 

15 

DA;, 

- 

♦ 1 , n 0 


i) A I 

- 

D A D - 1 . DO 


MAC = DMi+DAH ♦ DAI*0AL 

pn(; = pn IURPH ♦ DAI*D n L 

l)PB = OAC/0PC 

DAL = DAM 

i)«L « i)PH 

n Ah a I) AC 

DOH a ['PC 

I» ( nFA ,1 g.O.UO) on TU 20 

IFfnARSUDFP-DFAl/OFB). IE. 1.0-14) GO TO 20 
OF A * DFb 
GO TO 15 

20 YORMX * DPI*DFB*DEXP(-OX*OX/2.DO) 

if ( nx,LF,3,00) YURMX « 0,500 - YOPMX 
I F ( n Z » GT , 0 , DO ) VQRMX ■ 1.00 - YORHX 

RCTi'RM 

fcND 
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DMttt'LF PRFCJSION ^UNCTION OIOCH(OK) 

IMPLICIT REAl*e<D> 

DVmnx 

PTERM* 1 . DO 

IF (r»x) 1 * 1#2 

1 f)i r.r.Mso.DO 

RF ri'RM 

2 IF(PY-lft.DO) 3#3 #4 

3 mEPM«DTERH*OV 
DY« [)¥♦!'. DC 

GT TO 2 

4 DLOGH. ( DY-,500)* OLOG< DY > oDYM*. DO/ < 12,D0*DY)-1 ,00/ < 360, D0*DY4*3 ) 

1 +1,00/ ( 1260.00*DY**5) -1,00/ 1 16S0*D0*DVM7> ♦,91093853320467300 
2-r)LPG(DTERM) 

rfthru 

EMO 
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